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ABSTRACT 

This Letter analyzes 3-dimensional simulations of Kerr black hole magnetospheres 
that obey the general relativistic equations of perfect magnetohydrodynamics (MHD). 
Particular emphasis is on the event horizon magnetosphere (EHM) which is defined 
as the the large scale poloidal magnetic flux that threads the event horizon of a black 
hole (This is distinct from the poloidal magnetic flux that threads the equatorial plane 
of the ergosphere, which forms the ergospheric disk magnetosphere). Standard MHD 
theoretical treatments of Poynting jets in the EHM are predicated on the assumption 
that the plasma comprising the boundaries of the EHM plays no role in producing the 
Poynting flux. The energy flux is electrodynamic in origin and it is essentially conserved 
from the horizon to infinity, this is known as the Blandford-Znajek (B-Z) mechanism. 
To the contrary, within the 3-D simulations, the lateral boundaries are strong pistons 
for MHD waves and actually inject prodigious quantities of Poynting flux into the 
EHM. At high black hole spin rates, strong sources of Poynting flux adjacent to the 
EHM from the ergospheric disk will actually diffuse to higher latitudes and swamp any 
putative B-Z effects. This is in contrast to lower spin rates, which are characterized 
by much lower output powers and modest amounts of Poynting flux are injected into 
the EHM from the accretion disk corona. 
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De Villiers eTail j|2003l. 



simul ations disc uss ed i n 

20051); Ifflrose et all (|2004T ): 
Krolik et all (|2005l ); lHawlev and Krolikl (|2006l ) (HK, 



hereafter) offer an important virtual laboratory for study- 
ing the physics of black hole driven jets. These simulations 
readily evolve to a configuration with a net accreted poloidal 
magnetic flux that is trapped within the accretion vortex 
or funnel. This region is the black hole magnetosphere and 
it supports a jet dominated by electromagnetic en ergy for 
rapidl y rotating black holes. In a previous paper, IPunslvl 
(2007), it was shown that the phys ics conducive to an 
ergosp heric disk (first described in IPunslv and Coronitil 
(1990)) existed in the high spin simulation, a/M = 0.99 
(where the black hole mass, M, and the angular momentum 
per unit mass, a, are in geometrized units), known as KDJ. 
The ergospheric disk jet (EDJ) is launched from the plasma 
near the equatorial plane of the ergosphere. The event 
horizon magnetosphere (EHM) is comprised of poloidal flux 
that threads the event horizon (see the middle frame of 
figure 1). This paper explores the sources of Poynting flux 
in the EHM in the context of the 3-D simulations. In the 
EHM, there are four possible sources of Poynting flux that 
follow from energy conservation. 

(i) Energy can be transferred from the plasma to the elec- 



tromagnetic field. However, the EHM is virtually evacuated 
of plasma due to the centrifugal barrier and there is little 
energy to transfer to the field in perfect MHD. 

(ii) Any other sources must be surface terms at the 
boundary of the EHM. The first possibility is the electrody- 
namic component that is associated with the event horizon 
boundary surface. This is customarily called the Blandford- 
Znajek (B-Z) effect iBlandford and Znaiekl (|l977l ). 

(iii) The second boundary surface, which is associated 
with the EDJ, is the equatorial accretion flow in the er- 
gosphere. 

(iv) Outside the ergosphere, the EHM is bounded by the 
accretion disk and the disk corona. A source on this bound- 
ary will be called a coronal piston (see figure 2). 

It has been typically assumed that a Poynting jet in 
the EHM is confirmation of a B-Z process at work. How- 
ever, it will be shown that the two simulations, KDJ with 
a/M = 0.99 and KDH with a/M = 0.95, indicate that this 
can be a misleading interpretation in gene ral. The energ y 
output of KDJ is dominated by the EDJ IPunslvl l|2007h . 
KDH provides an interesting contrast to KDJ, because the 
EDJ is not nearly as powerful as it is in KDJ in the late 
time data slices. In KDJ, the Poynting flux in the EHM at 
large distances from the hole, r ~ 100M, is dominated by 
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Figure 1. The Poynting flux in KDJ. The color bar is in code 
units. The rows are in chronological order, t= 9840 M (top), 
t=9920 M (middle) and t= 10000 M (bottom). There is no data 
clipping, saturated regions are white. The region r < r^ n = 
1.203M is black and r+ = 1.141M. 



the energy flux from the EDJ that is gradually diffusing to 
higher latitudes within the funnel as it propagates outward. 
For KDH, the coronal piston that was discussed in HK in 
the context of driving the enormous mass flux in the "fun- 
nel wall jet" also injects significant amounts of Poynting flux 
into the funnel as a second order effect. 



1 THE 3-D SIMULATIONS 

J. Krolik and J. Hawley have generously shared the data for 
the last three time slices of KDH and KDJ, at t = 9840 M, 
t= 9920 M and t=10000 M. The simulations are performed 
in the Kerr metric (that of a rotating, uncharged black hole) , 
g M „. Calculations are carried out in Boyer-Lindquist coordi- 
nates {r,6,(j),t). The reader should refer to HK and refer- 
ences therein for details of the simulations. We only give a 
brief overview of the most relevant details. The initial state 
is a torus of gas in equilibrium that is threaded by concen- 
tric loops of weak magnetic flux that foliate the surfaces 
of constant pressure. The magnetic loops are twisted az- 
imuthally by the differentially rotating gas. This creates sig- 
nificant magnetic stress that transfers angular momentum 
outward in the gas, initiating a strong inflow that is per- 
meated by magneto-rotational instabilities (MRI). The end 
result is that after t = a few hundred M, accreted poloidal 
magnetic flux gets trapped in the accretion vortex or funnel 
(with an opening angle of ~ 60° at the horizon tapering to 
~ 30° - 35° at r > 20M). This region is the black hole mag- 
netosphere and it supports Poynting flux. The surrounding 
accretion flow is very turbulent. The strong transients die 
off by t = 2000 M, so the late time data dump s are the most 
physically relevant lHawlev and Krolik! (2006). 

Numerically, the problem is formulated on a grid that is 
192 x 192 x 64, spanning r m < r < 120M, 8.1° < 9 < 171.9° 
and < <f> < 90°. The inner calculational boundary, n n , is 
located close to, but just outside of the event horizon, r+, 
where the coordinates are singular. The <f> boundary condi- 
tion is periodic and the 9 boundary conditions are reflective. 
Zero-gradient boundary conditions are employed on the ra- 
dial boundaries, where the contents of the active zones are 
copied into the neighboring ghost zones. MHD waves prop- 
agate slower than the speed of light, therefore the gravita- 
tional redshift creates a magneto-sonic critical surface out- 
side of r+ from w hich no MHD wave can traverse in the 
outward direction IPunslvl (|200ll ). The philosophy was to 
choose ri n to lie inside the magneto-sonic critical surface, 
thereby isolating it from the calculational grid. There are 
also steep gradients in the metric derived quantities as r+ 
is approached. This is handled by increasing the resolution 
of the grid near r\ n with a cosh distribution of radial nodes. 
The validity of the numerics of this meth od was verified, 
near r in , in iDe Villiers and Hawlevl ([2003) by comparing 
simulations to simple analytic solutions. Even so, the simu- 
lations are closely monitored to look for unnatural bound- 
ary reflections. We also note the 3-D simulations in Kerr- 
Sc hild coordina t es (w hich are nonsingular on the horizon) 
in iFragile et 3 (120071 ). To test the code, they ran simula- 
tions of magnetized tori that were initiated from identical 
input parameters to those used by Hawley et al. In the words 
of C. Fragile (private communication), the results were "re- 
markably similar." Even though this was only verified for 
a/M=0.9, it is compelling. Consequently, for the purposes 
of this study it was concluded that the numerics were reli- 
able inside the ergosphere. 



2 THE KDJ SIMULATION 

The conservation of global, redshifted, energy flux, de- 
fined in terms of the stress-energy tensor, is simply, 
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9(V=gT l ")/9(i 1 ') = (|Punslvll2007h . The four-momentum 
— has two components: one from the fluid, — (T t ")fl u id, 
and one from the electromagnetic field, — (T t !/ )EM. The quan- 
tity g = — (r 2 + a 2 cos 2 8) 2 sin 2 is the determinant of the 
metric. The integral form of the conservation law arises from 
the trivial integrat ion of the partial differential expression 
iThorne et a 1 (1986). It follows that the poloidal components 
of the redshifted Poynting flux are S e = —y/—g (T t e )-EM and 
5"" = —y/—g (T{)-em- We can use these simple expressions 
to understand the EHM Poynting jet in KDJ. Figure 1 is a 
plot of S r in KDJ viewed at three different levels of magni- 
fication, at t= 9840 M (top row), t= 9920 M (middle row) 
and the bottom frames are at t = 10000 M. Each frame is 
the average over azimuth of each time step. This greatly re- 
duces the fluctuations as the accretion vortex is a cauldron 
of strong MHD waves. The individual (f> — constant slices 
show the same dominant behavior, however it is embedded in 
large MHD fluctuations. The left hand column show s strong 
beams of S r coming from near the black hole. In IPunslvl 
(|2007l V it was shown that the source of these beams was 
S 6 that was created by the ergospheric disk. The base of 
the EDJ emerging from the dense equatorial plasma is well 
resolved in both r and 9 with w 20 and > 30 grid zones, 
respectively. 

In this paper, we turn our attention to the propagation 
of individual flares from the ergospheric disk out to the outer 
calculations! boundary at r = 120M. Even though the time 
sampling is very coarse in the data dumps (At = 80M), 
we can understand the propagation of the EDJ because 
of the wide angle views available in the right hand col- 
umn of figure 1. We track the EDJ evolution by identi- 
fying the strong knots or flares in figure 1 based on the 
following reasoning. As discussed in IPunslvl (|200lf ). the S r 
flares will propagate at the speed of an MHD discontinuity 
as modified by the plasma bulk flow velocity. The plasma 
near the edge of the vortex has accelerated to v r > 0.9c 
by r= 30 M. So the flares of S r should propagate radially 
at Vfiare ~ c for r > 30M. Without having the benefit of 
the detailed time evolution, this upper bound is the best 
estimate that we can make for Vfi are . First, consider the 
strong knot, "C," at t=10000 M in the right hand frame. 
Label the outer radial extent of knot C at t = 10000 M 
by r +c (t = 10000A/) = 100.9A/ and inner radial edge by 
r- C (t = 10000M) = 65.9M. Translating this MHD discon- 
tinuity back in time to t = 9920 M is equivalent to a radial 
displacement V flare At = V> Jore (-80M/c) > - 80M. Thus 
at t= 9920M, knot "C" should extend from the ergospheric 
disk to r +c (t = 9920M) > 20.9M. This is verified by the red 
patches in the middle frame at t=9920 M. 

Next consider the strong knot, "A," at t= 9840 M in 
the right hand frame. Label the outer radial extent of knot 
A at t = 9840 M by, r +A {t = 9840M) = 58.9M and in- 
ner radial edge by r_ A (t = 9840M) = 22.0M. Time trans- 
lating this feature to t = 9920 M implies that r-A{t = 
9920A4") < 102. 0M, so it must be visible near the edge of 
the right hand frame at t = 9920 M. Furthermore, unless 
the flare is propagating inordinately slowly, Vfi are < 0.75c, 
r+A(t — 9920M) will be beyond the outer boundary of the 
plot. There is only one plausible feature at t =9920 M. Fi- 
nally, there is a strong flare "B" that is emerging from the 
ergospheric disk at t = 9840 M in the left and middle frames. 
Thus, at t = 9920 M, some portion of the flare must be 



within 80 M of the black hole, hence the identification of 
"B" in the right hand frame. 

Figure 1 demonstrates a type (iii) source in the notation 
of the introduction. At r w 1.5M - 2.0M, the EDJ enters 
the EHM from the periphery (left hand frames). The EDJ 
gets quickly linked into the EHM because the ergopsheric 
disk magnetosphere in KDJ is comprised of small patches 
of twisted vertical flux that become intertwined with the 
large -scale flux in the EHM on scales of ~ 1M - 2M (|Punslvl 
2007). After this rapid injection, the S r in the EDJ keeps 
spreading towards the pole as it propagates outward. At 
the time steps that were made available to this author, 
the EDJ is the predominant source of 5"" in the EHM. By 
r » 100M, the EDJ is flooding the EHM, even close to the 
polar axis (it should be noted that the total energy flux is 
larger than S r , « 20% is in mechanical form in the EHM 
at r ~ 100M). The relevance to this discussion is that the 
EHM is inundated with 5"" that was not created on field 
lines that thread the horizon, but on flux entrapped within 
the equatorial accreting plasma. The slow diffusion of S r 
poleward at r > 30M is most likely regulated by numer- 
ical diffusion. This might seem like a problem from a nu- 
merical point of view, but physically this is not nearly as 
much of a concern from a qualitative standpoint. Perfect 
MHD is just a simple tractable method of dealing with the 
plasma physics. A realistic, high temperature, jet plasma is 
likely to have anomalous res i stivity from a varie ty of sources, 
ISomov and Oreshinal (l200(il ) ; iTreumannl (| 200ll). and the dif- 
fusion of field energy should naturally occur. The simula- 
tion cannot accurately describe the diffusion rate. However, 
qualitatively speaking it indicates that if the jet propagates 
extremely far from the hole (r 2> 120A/), regardless of the 
exact details of the diffusion microphysics, the EDJ energy 
flux is likely to get smeared out towards the polar region. 



3 THE KDH SIMULATION 

There is a relatively weak EDJ in KDH that is noticeable at 
t = 9840 M, but it is otherwise negligible. This circumstance 
allows for the detection of weaker sources of S r that would 
otherwise be swamped by a strong EDJ. Figure 2 is a plot 
of ST for these three time slices in chronological order, going 
from left to right. Each frame is the average over azimuth 
of the time step. The contours of the radial momentum flux 
due to mass motion, P r = ^/—gpU r U r (where p is the proper 
density and is the four velocity) , is overlayed in white in 
order to define the location of the "funnel wall jet," as was 
done in HK. The funnel wall jet is a shear layer between the 
accretion disk corona and the Poynting jet. It is a collimated 
sub-relativistic flow that transports most of the mass outflow 
in the jetted system. In HK, it was shown to be driven by 
the total pressure (gas plus magnetic) gradient in the corona 
that is oblique to the funnel wall boundary. The gas in this 
region is constrained from being pushed into the funnel by 
the centrifugal barrier. The component of pressure gradient 
that is parallel to the centrifugal barrier forces the flow to 
be squeezed outward as a shear layer. 

In figure 2, there is an almost one to one correspon- 
dence, between locations where P r of the funnel wall jet 
increases and sites where S r increases at the funnel wall 
boundary in the EHM. It is the coronal pressure and not 
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Figure 2. The coronal MHD piston is illustrated by these plots of the Poynting flux in KDH. The color bar is in code units. The frames 
are in chronological order, t= 9840 M (left), t=9920 M (center) and t= 10000 M (right). There is no data clipping, saturated regions 
are white. The overlayed white contours represent, P r , the radial momentum flux described in the text. At each coronal piston location 
there is a large pressure flare, see figure 3. Inside of r^„ = 1.403M (r_|_ = 1.312A/) is black. 




Figure 3. A closcup of the coronal piston at t =10000 M is 
depicted in these plots of the total pressure density (gas plus 
magnetic), expressed in terms of the gas pressure P g and the 
Faraday tensor as \/—g{P g + F^ F^v/lQir). The left frame is 
averaged over rfr and the right frame is at = 49.2°. The color bar 
is the strength of the total pressure in code units. High pressure 
regions of the corona are saturated and appear white. The force 
associated with the pressure gradient is indicated by the white 
arrows. The overlayed white contours represent, P r , the radial 
momentum flux as in figure 2 and the black contours represent 
S r . 



the outgoing Poynting jet (OPJ, hereafter) that drives the 
funnel wall jet. Figure 3 shows a strong flare in the total 
pressure (gas plus magnetic) at t =10000 M. The pressure 
gradient seems to provide the accelerating force that drives 
P r . The flare appears to be a high pressure loop emerging 
from the corona, as evidenced by the right hand frame of fig- 
ure 3. The loop location and topology are inconsistent with 
the high pressure feature being injected into the corona from 
the funnel interior. The magnetic pressure in the corona ac- 
tually exceeds the magnetic pressure in the EHM at these 
intermediate radii, lOAf < r < 30M, precisely the region 
where the flares tend to occur. Table 4 of HK indicates that 
S~ = f r=r S r d6d<t>dt = J Sl n d6 = (2.79/4.26)P jet , where 
Pjet is the total energy transported tor = 120Af by the fun- 
nel wall jet. Table 4 also indicates that J r _ 120M S r d9d<f>dt 
— (1.46/2.79) Sin- Therefore, by conservation of energy, a 



putative B-Z effect can provide at most (2.79— 1.46)/4.26 = 
0.31 of Pj e t, and is too feeble to drive the funnel wall jet. 

Most of Sin is directed into the equatorial accretion 
flow and never reaches the plasma in the OPJ. To see this, 
we summarize figure 5 of HK. There are broad relative 
maxima in S|" n at 9 < 41.3° (9 > 137°), but Sl n decreases 
rapidly with 9 and S- n < when 41.3° < 9 < 54.5° 
(120.4° < 9 < 137°). The two strongest maxima in S r „ are 
in the equatorial accretion flow at 9 w 80° (9 « 100°) and 
this energy never reaches the OPJ. Therefore, it is unclear 
how the region 41.3° < 9 < 137° is related to the B-Z mech- 
anism. In order to understand the global energy budget, we 
want to know how much Sin reaches the OPJ. Based on the 
available data, we ca n generate an uppe r bound. First, we 
define the OPJ as in Hiro se et a 3 (|2004l ) by the magnetic 
dominance condition, defined in terms of the Faraday field 
strength tensor as F^" F nv / '(lGirph) > 1 (ft is the enthalpy 
per unit mass) and U r > 0. Virtually all the S r that reaches 
r — 120M in the available KDH time slices flows in the OPJ 
as defined above. Inspecting our three time slices, the OPJ 
initiates at r m i„ < 3M, with an angular extreme near the 
base given by 9 ma x ~ 55° {9 min » 125° ) . in agreement with 
the ti me averaged data in figure 11 of lHawlev and Krolikl 
(2006). The OPJ is collimated, i.e., d9 max /dr < for 
9 < 90° (d9 min /dr > for 9 > 90°). In each of the 
time slices, we can compare J S r d9d<j) near ri„ above the 
equatorial accretion flow to J S r d9d(f> in the OPJ to assess 
how much Si n reaches the OPJ. The simulations have 
significant fluctuations. Averaging over r smooths out the 
spatial fluctuations and is a more reliable diagnostic than 
computing fluxes at a single radius. Unfortunately, there 
is not enough data to perform a meaningful time average. 
In each time slice we find 5 g 7M J^ M dr Jg_^° 10 S r d9d(f> > 

iu C *■ i:::rf *s££" sr ded <t> > 

337 Ism dr le-e ^ d8d(j). The 6 inequalities tend to 
indicate that Sin at 9 < 70° (9 > 110°) is larger (the 
average excess is 20%) than the total S r flowing through 
the OPJ at r w 7M, 9 < 50° (6 > 130°) and some of the 
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Sin at 9 < 70° is absorbed by the accretion disk and corona 
at r 5? 7M, 9 > 50° (9 < 130°). This excess is apparently 
large enough to persist in spite of significant temporal 
variations. Conversely, the probability of the hypothesis 
that random temporal fluctuations have conspired to mask 
what is actually a conserved energy flow from ri n , 9 < 70° 
(6> > 110°) to the OPJ atrsi 7M is rejected at the (0.5 



per hemisphere) 1 — (0.5) 



0.9844 significance level. 



Integrating the area under the plot of 57„ in figure 5 of 
HK, one finds that only 39% of the total Si n originates in 
the angular range 8.1° < 9 < 70°, 110° < 9 < 171.9°. Thus, 
the data tends to show that Q.39S in is a loose upper limit 
to the maximum amount of Si n that can reach the OPJ. 
Applying this to the results in table 4 of HK implies that 
> [1.46 - (0.39)(2. 79)]/1.46 = 0.26 of the 5"' in the OPJ at 
r — 120M is created at r > ri„ during the course of the sim- 
ulation. Furthermore, integrating over the funnel cross sec- 
tion at both ri n and across the flares, [Jjj~ ™ a S r d9d(j>] flare > 
1.5(1/0.597M) J* M dr J^™1 S r d9d<f> in the last two time 
steps (there is a weak EDJ at t =9820 M that skews our 
results, so it is omitted from this analysis). If one includes 
the mechanical energy flux as well, there is consistently 
more than twice the total energy flux in the strong flares 
then there is total energy flux through Vi n , 9 < 70°. This 
suggests that coronal injection sites are sufficiently strong 
to make up this apparent > 26% deficit in S T . 

One might be concerned that there are only « 8 -10 an- 
gular zones between the coronal piston and the funnel and 
this leads to significant numerical diffusion. From a numer- 
ics point of view this is much more of a concern than from a 
physical point of view. The MHD code is just a simple ap- 
proximation to any real turbulent plasma state. The turbu- 
lent corona is likely to have an anomalous resisti v ity and dif- 
fusion should occur ISomov and Oreshinal l|2000r i: iTreumannl 
l|200lf ). The rate of diffusion cannot be determined by this 
simulation. However, the qualitative idea that a strong flare 
in coronal energy can in principle reach the funnel interior 
is strongly indicated. 



4 DISCUSSION 

The philosophy of this paper is not that the simulations 
of idealized magnetized tori gives us a direct picture that 
can be applied to AGN central engines. They are treated 
only as virtual laboratories to see what effects might self- 
consistently occur near a magnetized black hole. These sim- 
ulations have led us to a new realization, the boundaries of 
the EHM should be dynamic and are not likely to be passive 
boundary surfaces for the magnetic field. It was shown that 
electrodynamic energy flux can arise in the EHM as a result 
of sources radiating energy from the lateral boundaries. Even 
if the EHM can be construed as "force- free," the dynamics 
of the lateral boundaries are determined by strong inertial 
forces that should make them strong MHD pistons. This 
circumstance was not anticipated in theoretical treatments 
of ele ct rodynam i c jets in the EHM iBlandford and Znaiekl 
1 19771 ) ; |Phinnevl (| 19831 ). The fact that electromagnetic en- 
ergy can come into the EHM from the side goes right to the 
heart of the assumptions in the B-Z solution. The B-Z solu- 
tion is the perfect MHD solution in which energy conserva- 



tion reduces to Poynting flux conservation from the horizo n 
to a relativistic wind at asymptotic infinity [Phlnnevl l|l983h . 
From this condition, the parameters of the field are uniquely 
determined for a given poloidal field distribution, in partic- 
ular the field line angular velocity, Qf, and the total elec- 
tromagnetic energy output from the black hole, J S r d9d<f>. 
If there are strong sources of Poynting flux along the lat- 
eral walls of the EHM, the spacetime near the event horizon 
can not adjust the system to enforce the B-Z field param- 
eters within the EHM. This is a direct consequence of the 
fact that the plasma near the event horizon in the EHM can 
not effectively react back on the outgoing wind or jet and 
modify its electromagnetic properties becaus e of the gravita - 



tional redshifting of the MHD characteristics IPunslvi (|200lT ) ; 
iPunslv and BinJ i|2004T ). The plasma near the horizon in the 
EHM will passively accept any field para meters impose d by 
the EDJ and the accretion disk corona (IPunslvi 1200 ll ). As 
such, in a general astrophysical context, the basic parame- 
ters such as H,f and J S r d9d(f> are indeterminant. Of course, 
this does not preclude the possibility of an MHD numerical 
system evolving tow ards B-Z, if the num erical problem is 
properly constructed iKomissarov I (|2004f) . From a physical 
point of view, the strong forces that are responsible for com- 
pressing the flux down to the horizon still reside in the " fun- 
nel walls," rendering the lateral boundaries as strong MHD 
pistons. In a realistic astrophysical setting, inertial forces in 
the lateral boundaries are likely to play an important role, or 
even a dominant role, in the determination of the jet power 
from the EHM. 
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